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Top quark phenomenology is one of the main fields of investigation in theoretical and experi- 
mental particle physics. The experiments at the Tevatron accelerator and, in future, at the LHC will 
allow one to perform improved measurements of the top properties, such as its mass, thanks to the 
production of large amounts of tt pairs. 

In this paper we study bottom quark fragmentation in top quark decay (t — > bW), which is 
responsible of one of the largest contribution to the uncertainty on the top mass measurement at 
the Tevatron [jl], ^ and the LHC [|3p. In particular, in the analysis of Ref. [Q] the top quark mass is 
determined using at the LHC final states with leptons and //y's, where the leptons come from the 
W decay W — > Iv, and the J/y's from the decay of a ^-flavoured hadron B. In [Q] the PYTHIA 
event generator [Q| was exploited, and the error on m t was estimated to be Am, ~ 1 GeV, with 
^-fragmentation being the largest source of uncertainty. In [|||, the invariant mass mst, yielded by 
the HERWIG [f7|] event generator, was used to fit m t , and the impact of matrix-element corrections 
to the simulation of top decay [||] was investigated. 

Bottom quark fragmentation in top decay was studied in []|, [To| ], following the method of 
perturbative fragmentation functions [JTTJ]. The NLO b-quark energy spectrum is expressed as the 
convolution of a coefficient function, describing the emission of a massless parton, and a pertur- 
bative fragmentation function D(m,b,iip), associated with the transition of a massless parton into 
a massive b. D(mb,jj,f) follows the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) evo- 
lution equations [|l3|, [14]], which can be solved once an initial condition at a scale /Iof is given. 
The initial condition of the perturbative fragmentation function, first computed in [pT|], was proved 
to be process-independent in Jl2|]. Solving the DGLAP evolution equations we can resum the 
large \n(nv} /mfy which appears in the NLO massive ^-spectrum (collinear resummation). Both 
the top-decay coefficient function, computed in |9p, and the initial condition D(mb,Hof) present 
terms which become large when the ft-quark energy fraction x/, approaches 1 , which corresponds 
to soft-gluon radiation. Soft contributions in the initial condition (process independent) and in the 



coefficient function (process dependent) were resummed in the NLL approximation in Jl2| ] and 



[10], respectively. In order to predict the spectrum of ^-flavoured hadrons, perturbative calcula- 



tions need to be supplement by non-perturbative models. In []9[ [l(J, the B-hadron spectrum in top 



decay was presented, after fitting a few hadronization models to SLD [15] and ALEPH [[16] data. 



Following the lines of Q17[], in this paper we would like to address ^-fragmentation in top decay, 



using the PYTHIA and HERWIG event generators. As discussed in Q17J1, PYTHIA and HERWIG 
simulate multiple radiation in top decay in the soft or collinear approximation, and are provided 
with matrix-element corrections [|8, 18] to allow hard and large-angle emission. The hadronization 



mechanism is simulated by the string model [|19p in PYTHIA, and by the cluster model Q20[ ] in 
HERWIG. 

For the sake of a reliable prediction of the B-energy distribution in t — > bW, we need to use 
models and parametrizations which are able to describe well the B-spectrum at e + e~ machines. We 



consider ALEPH Q16[], OPAL Q2 1[ ] and SLD []15J data on the B energy fraction xb in Z — > bb events, 



where x# is defined as follows: 



2-Pb-Pz m 
xb = ~i — , (1) 



with pz and ps being the Z and B momenta, respectively. 
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Table 1: Parameters of HERWIG and PYTHIA hadronization models tuned to e + e data, along with the 
X 2 per degree of freedom. 



HERWIG 


PYTHIA 


CLSMR(l) = 0.4 




CLSMR(2) = 0.3 


PARJ(41) = 0.85 


DECWT = 0.7 


PARJ(42)= 1.03 


CLPOW = 2.1 


PARJ(46) = 0.85 


PSPLT(2) = 0.33 




% 2 /dof =222.4/61 


X 2 /dof = 45.7/61 



We use the versions HERWIG 6.506 and PYTHIA 6.220, and find that the default parametriza- 
tions are unable to fit such data, as one gets # 2 /dof = 739.4/61 (HERWIG) and £ 2 /dof = 467.9/61 



(PYTHIA). As in [ 17], we tune the cluster and string models to the data, while we leave unchanged 
the parameters of HERWIG and PYTHIA which are related to the perturbative phase of the parton 
showers. Our best fits are summarized in Table [I]: for PYTHIA we are able to find a parametriza- 
tion which is able to reproduce well the data (# 2 /dof = 45.7/61); HERWIG, even after the fit, is 
still not able to describe the jcg-spectrum very well, but the comparison is anyway much better than 
with the default parameters (# 2 /dof = 222.4/61). We have also checked that the parametrizations 
in Table [l] work well for the new model implemented in PYTHIA 6.3 The HERWIG and 

PYTHIA spectra, before and after the fit, along with the experimental data, are presented in Fig- 
ures [I] and|2[ For the sake of comparison, we also show the Ag-spectrum yielded by the NLO+NLL 
calculation of Ref. [|12|], convoluted with the Kartvelishvili hadronization model [23]: 

D n P(x; 7 ) = (l + 7 )(2 + 7)(l-x)x^ (2) 

We fit the model (^) to the data in the range 0.18 < xb < 0.94, to avoid the regions at small and 
large xb, where, as pointed out in [JTJ], the resummed calculation yields a negative distribution. 
Setting mz = 91. 118 GeV, nit = 5 GeV and A™. = 200 GeV in the perturbative calculation, we 
get 7 = 17.178 ± 0.303 From Figs. [I] and | we learn that default HERWIG and PYTHIA are far 
from the data. After the tuning, PYTHIA reproduces the data quite well, while HERWIG yields a 
broader distribution; the resummed calculation is consistent with the data, and # 2 /dof = 46.2/53 
from the fit. 

Using the parametrizations in Table [I], we can predict the B-energy distribution in t —* bW, 
which will be expressed in terms of the variable 

xb = z — , (3) 

1 — w mf 

where p t is the top momentum and 1/(1 —w) is a normalization factor, with w = 1 — m^/m 2 + 
m 2 h /m 2 . In Fig. |3| we present the B-spectrum in top decay according to HERWIG, PYTHIA and 



the resummed calculation of M, 10], convoluted with the Kartvelishvili model. The comparison 



exhibited in Fig. || shows similar features to Figures [j] and ^, and reflects the quality of the fits to 
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Figure 1: Data from LEP and SLD experiments, compared with the NLO+NLL calculation convoluted with 
the Kartvelishvili model (solid) and HERWIG 6.506, using the default parametrization (dashed) and our 
tuning (dotted). 




0.2 0.4 0.6 0.8 1 



X B 

Figure 2: As in Fig. [IJ but comparing data and the NLO+NLL calculation with default (dashed) and tuned 
(dotted) PYTHIA 6.220. 



the e + e~ data: PYTHIA reproduces the peak of the resummed calculation rather well, while it is 
below the NLL prediction at xb < 0.7, and above at x# > 0.9. HERWIG is below the resummed 
spectrum in most the x^-range, and above it only at large xg. 

Finally, we wish to present results in Mellin moment space, making use of the data reported by 
the DELPHI Collaboration in [24] on the moments of the B cross section in e + e~ processes. From 
the point of view of resummed calculations, working in moment space presents several advantages: 



4 



B-hadron production in top quark decay 



Gennaro Corcella 
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X B 

Figure 3: B-hadron spectra in top decay, for m t = 175 GeV, according to a NLO+NLL computation convo- 
luted with the Kartvelishvili model (solid line), HERWIG (dashed) and PYTHIA (dotted). 



in Af-space, convolutions become ordinary products and the relation between parton- and hadron- 
level moments is just o$ = <7^D^ P , where D^ p is the non-perturbative fragmentation function in 
Mellin space. Therefore, there is no need to introduce a functional form for the hadronization 
model in .xg-space. Also, N- spectra are well defined, and do not present the problems of the 
xg-results, which are negative at small or large x B . 

In Table || we compare the DELPHI moments with the ones given by the tuned versions of 
HERWIG and PYTHIA, and the predictions for top decay. We also quote the results yielded by the 
NLL calculations of Refs. []9[ [T(J, extracting D^ p from the data. As for e + e~ — ► bb processes, the 
moments given by HERWIG and PYTHIA are consistent, within the error ranges, with the DELPHI 
ones. It is interesting that HERWIG is compatible with the DELPHI moments, even though it was 
only marginally consistent with LEP and SLD data in xg-space. The results for top decay have 
similar features to the xg-spectra: PYTHIA is very close to the NLL calculation, which uses D^ p 
extracted from the DELPHI data, while HERWIG, whose predictions are shifted toward larger xg, 
gives larger moments. 

In summary, we reviewed recent results on ^-flavoured hadron production in top quark decay. 
We tuned HERWIG and PYTHIA to LEP and SLD data and presented results on the B-hadron 
spectrum in top decay in both x B and moment spaces. In fact, fitting the cluster and string model 
turned out to be necessary to reproduce the e + e~ data. The results were also compared with re- 
summed calculations, based on the fragmentation function formalism. We think that our analysis 
and fits can be useful to improve the present understanding of &-quark fragmentation in top quark 
decay. It will be very interesting to investigate how the tuned versions of HERWIG and PYTHIA 
fare with respect to other observables. For example, it may be worthwhile reconsidering the studies 
in Refs. [|], ||] with the parametrizations which we have proposed, and estimate the contribution of 
^-fragmentation to the uncertainty on the top quark mass reconstruction. 
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Table 2: Moments a§ from DELPHI [24|, and moments in e + e annihilation and top (f) decay, using NLL 
resummed calculations, HERWIG (HW) and PYTHIA (PY). 





(x) 


(x 2 ) 




i A \ 

(x A ) 


e + e data 


0.7153±0.0052 


0.5401 ±0.0064 


0.4236±0.0065 


0.3406±0.0064 


e+e- NLL a* 


0.7801 


0.6436 


0.5479 


0.4755 


U N 


0.9169 


0.8392 


0.7731 


0.7163 


e+e- HW o$ 


0.7113 


0.5354 


0.4181 


0.3353 


e+e- PY o# 


0.7162 


0.5412 


0.4237 


0.3400 


f-dec. NLL T b N 


0.7883 


0.6615 


0.5735 


0.5071 


f-dec. NLL T% = r^D^ 


0.7228 


0.5551 


0.4434 


0.3632 


f-dec. HW 


0.7325 


0.5703 


0.4606 


0.3814 


f-dec. PY T% 


0.7225 


0.5588 


0.4486 


0.3688 
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